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摘要 : 数据 融合 的 目的 是 合并 或 者 综合 多 个 数据 集 ， 成 为 一 个 完整 的 性 能 更 好 的 数据 
集 。 基 于 广义 延 拓 插值 外 推 法 提出 建立 两 组 或 两 组 以 上 方程 解 集 和 数据 序列 融合 优化 的 数学 模 
型 ， 称 之 为 广义 融合 法 ， 此 方法 能 充分 发 挥 数据 和 量 测 资 源 的 潜力 ， 运 用 创新 的 集成 和 互 校 算 
法 ， 解 决 数据 融合 面临 的 难题 。 介 绍 了 广义 融合 数学 模型 的 建立 和 求解 方法 的 主要 步骤 ; 对 模 
型 中 权 值 的 分 配 原则 做 了 分 析 ; 最 后 应 用 于 卫星 导航 数据 处 理 领 域 。 实 验证 明 该 方法 实用 性 
好 ， 数 据 平稳 性 好 ， 外 推 稳定 度 高 ， 有 望 成 为 一 种 新 的 解决 数据 融合 问题 的 普 适 求解 方法 。 
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T 数据 融合 (也 叫 信息 融 合 ) 的 概念 始 于 20 世纪 70 年 代 ， 至 80 年 代 其 技术 取得 了 一 定 的 进步 和 发 
` 展 ， 到 90 年 代 不 少数 据 融合 技术 的 研究 成 果 已 得 到 实际 应 用 ， 取 得 了 理想 的 效果 '' 1。 数据 融合 的 主 
ME 要 原理 是 用 多 个 传感器 ( 同 质 或 异 质 的 ) 对 同一 对 象 进行 描述 ， 得 到 该 对 象 的 多 源 观 测 数 据 ， 将 这 些 
m 数据 进行 集成 ， 以 获得 对 被 控 对 象 的 更 好 的 理解 号 。 
~ 多 传感器 信息 融合 技术 是 近 几 年 发 展 起 来 的 一 门 实践 应 用 技术 ， 是 针对 一 个 系统 使 用 多 种 传感器 
这 一 特定 问题 而 展开 的 一 种 关于 数据 处 理 的 技术 5 。 对 于 多 传感器 系统 ， 由 于 信息 的 多 样 化 和 复杂 
性 ， 因 此 ， 信 息 融 合 方法 的 基本 要 求 是 具有 强大 的 处 理 能 力 ， 以 及 较 高 的 计算 速度 和 精度 。 人 工 智能 
的 发 展 使 信息 融合 技术 走向 智能 化 、 集 成 化 。 

数据 融合 的 目的 是 合并 或 合成 多 个 数据 集 ， 成 为 一 个 完整 的 性 能 更 好 的 数据 集 。 数 据 融合 技术 的 
实现 ， 关 键 在 于 融合 算法 !#5 。 在 融合 求解 过 程 中 ， 不 是 仅仅 把 数据 或 解 集 简单 地 直接 组 合 ， 这 样 做 
效果 不 会 很 明显 ， 甚 至 最 后 得 到 的 结果 更 差 。 为 了 能 够 得 到 良好 的 融合 效果 ， 应 该 分 析 各 组 方程 解 或 
PT 多 组 数据 集 序列 的 特征 。 比 如 一 些 数据 序列 的 随机 误差 较 小 但 偏差 较 大 ; 一 些 数据 的 随机 误差 较 大 但 
偏差 较 小 "1 。 这 样 就 可 以 将 随机 误差 较 小 的 数据 序列 的 相对 变化 量 通过 一 定 的 方法 融合 到 随机 误差 
较 大 的 数据 序列 中 ， 能 获得 随机 误差 较 小 且 偏 差 也 较 小 的 数据 序列 。 本 文 基于 广义 延 拓 插 值 外 推 
法 中 提出 了 一 种 求解 一 类 特定 的 关联 群 方程 组 的 多 个 离散 数据 集 或 最 优 融 合 估 计 集 的 方法 ， 称 之 为 
广义 融合 方法 。 经 过 实验 证 明 ， 该 方法 实用 性 好 ， 数 据 平稳 性 好 ， 外 推荐 定 度 高 ,已 应 用 于 卫星 导航 
定位 数据 处 理 领 域 。 


1 广义 延 拓 插值 外 推 法 


广义 融合 是 基于 广义 延 拓 插值 外 推 法 提出 的 ， 因 此 ， 先 介绍 广义 延 拓 插值 外 推 法 。 广 义 延 拓 插 值 
外 推 法 的 基本 设计 理念 是 ” : 分 段 实 施 最 小 二 乘 逼近 ,在 每 段 的 端点 处 采用 插值 处 理 ， 在 段 内 采用 
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拟 合 处 理 ， 段 内 拟 合 通 近 处 理 时 ， 用 邻近 延 拓 域 
里 的 数据 进行 辅助 拟 合 处 理 ， 从 而 使 盘 近 多 项 式 


在 定义 域内 的 和 近 精度 升 高 ， 当 分 段 连接 后 , 在 
连接 端点 处 过 渡 平滑 ， 在 整 域 上 逼近 精度 比较 e OOT 


好 。 广 义 延 拓 最 小 二 乘 盘 近 方法 的 外 推 模型 兼 融 
了 插值 和 拟 合 两 大 功能 ， 具 有 极 好 的 逼近 精度 。 | 
外 推 示意 图 见 图 1。 | 

假设 已 知 历 元 时 刻 的 最 优 估计 值 六 入 | 
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之 前 的 数据 值 (各, 各，…, ia), 求 下 一 时 刻 。 0 mmm 
ta BIE ww， 这 是 一 个 典型 的 数据 外 推 问题 ， 图 1 广义 延 拓 外 推 示意 图 
可 建立 如 下 广义 延 拓 插值 儿 PHERI E : Fig.1 Generalized extrapolation diagram 
min/(a,, a,, a3) = Y [a, + at; + azt; -4]^ 
i-m 9 (1) 
s.t. d, t ast, 十 at, = X, 
p^ Jr, min/(a, , a; , a; ) 为 通 近 残 差 极 小 化 目标 函数 ，o a, a, 为 待 求 的 系数 ，s.t. 示 意 约束 条 件 。 
e 求解 模型 (1) ， 可 以 优化 获得 广义 延 拓 通 近 多 项 式 的 待 求 系数 w as. a4; 代 和 人 (2) 式 : 
x(t) = a, + at; + azt; (i=nt+l1,n+2, …), (2) 
c fi T DAPMRORÁRERERI c, Do TOME S, 
T JA LT VIA YY H- ` Yas 
一 2 广义 融合 方法 的 数学 模型 及 其 求解 方法 
~N 
O 2.1 广义 融合 方法 数学 模型 的 建立 


广义 融合 求解 模型 借鉴 了 广义 延 拓 插值 外 推 公 式 ， 即 在 i 历 元 时 ， 选 用 多 项 式 atat, tau, 作为 
逼近 函数 。 通 近 时 选用 一 段 历 元 内 的 数据 ， 其 中 把 已 获得 的 先 验 最 优 佑 计 值 (起 始 时 由 于 没有 先 验 最 
优 估计 值 ， 可 以 先 用 测量 值 代替 ) 作 拟 合 逼 近 处 理 ; 选用 当前 待 求 历 元 的 广义 状态 量 的 预 估 值 与 测量 
值 的 最 优 加 权 组 合 台 近 形式 作为 插值 约束 锁定 点 ; 以 逼近 多 项 式 系 数 w as. a, RUBUS EL, E, 等 
为 优化 求解 变量 ;以 允 近 函数 多 项 式 系数 a; a, a, IRÉ k, k, 以 及 广义 状态 量 等 待 求 量 的 区 
间 数 作为 变量 取 值 约束 条 件 。 这 些 约束 条 件 把 解 限制 在 一 定 的 范围 之 内 ， 从 而 提高 解 的 精度 。 

构造 广义 融合 最 优化 求解 模型 ， 求 a a 、as ， 使 满足 : 


n-2 
min/,(a,, 45, az) = Y [a + at; + azt; e. 
i-n-m 
k 2 ^ i k 2 2 
十 iLa 十 at, + azt, =w | T ola; * at, + azt, =x] (3) 
3 


> 
s.t. a; ta, , tau, Tx, as; k tk; =l 


n-l 


a Câ, a Ca, a Ea 

其 中 , Ca, a, 24) 为 待 优化 的 目标 函数 ， 优 化 的 目标 函数 既 可 以 是 线性 关系 也 可 以 是 非 线性 关 
系 ， 若 是 非 线性 关系 ， 则 广义 融合 方法 直接 采用 非 线 性 求解 算法 求解 ; *; 为 状态 量 的 测量 值 ; zx 
为 状态 量 的 最 优 估计 值 ; * 为 由 低 阶 状态 量 递 推 获得 的 状态 量 的 预 估 值 ; ;为 历 元 序号 ; k, k, 为 权 
系数 ,它们 的 取 值 大 小 影响 它们 各 自 的 残 差 对 目标 函数 的 贡献， 权 值 的 分 配 在 本 文 后 面 作 详 细 的 
介绍 ; x, 为 t 历 元 时 状态 量 的 测量 值 ; x, 为 1 历 元 时 状态 量 的 递 推 预 估 值 ，%, 为 4， 
历 元 时 状态 量 的 最 优 估 计 值 a,、 纪 了 、a 为 可 行 区 间 变 量 a, a5, a, 的 约束 区 间 数 ，41=[al、ai]， 
a =la 3]，@3=[a3、@3]。 广 义 融 合 模型 中 可 以 包含 非 线 性 状态 方程 、 非 线性 约束 方程 以 及 其 它 非 
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线性 测量 方程 ， 这 能 够 扩大 应 用 范围 ， 增 加 方程 的 元 余 性 和 连续 性 ， 提 高 求解 精度 。 
2.2 广义 融合 方法 的 求解 步骤 
2.2.1 求解 状态 量 值 
建立 残 差 极 小 化 优化 模型 并 加 约束 方程 ， 如 下 式 : 
minf(x,) = Y, (FG) = Y)? 
s.t. f(x,)-0 | 
HP, Fa) 为 状态 量 关系 式 ， 多 数 情 况 下 为 非 线性 函数 关系 ， 可 以 采用 非 线性 算法 直接 求解 ; x,; 为 
函数 状态 变量 ; Y, 为 测量 量 ; i 为 历 元 序号 ; n 为 当前 历 元 数 ; f(x, ) 为 约束 函数 关系 式 。(4) 式 可 以 
采用 非 线 性 算法 求解 ， 如 单纯 形 法 等 。 通 过 搜索 、 比 较 等 步骤 ,使 目标 函数 值 不 断交 近 最 优点 。 求 解 
方法 在 优化 过 程 中 不 受 方程 个 数 的 限制 ， 并 且 具 有 简单 直观 、 适 应 性 强 、 应 用 范围 广 、 解 决 速度 快 等 
寺 点 。 
2.2.2 求解 状态 量 的 预 估 值 c, 
利用 外 推 得 到 的 状态 方程 的 高 阶 状 态 量 可 以 提高 解 的 精度 和 相关 性 ， 低 阶 状态 量 可 以 是 与 状态 
相关 的 导数 值 或 微分 值 。 当 求解 i 历 元 时 刻 状态 量 的 预 估 值 *, 时 ， 由 于 已 经 得 到 1 历 元 时 刻 状态 量 
的 最 优 估计 值 x,_,( 在 初始 阶段 的 融合 ， 初 始 状 态 最 优 估 计 值 是 未 知 的 ， 可 用 前 几 个 时 刻 状态 的 测量 
ERE), ， 可 以 通过 公式 : 


(4) 


el Ka 


~oa : I: 
X, = Nil + x, At 二 PR deem , (5) 


递 推 获得 六 历 元 时 刻 状态 量 的 预 估计 值 x,。(5) 式 中 %, 为 4 历 元 时 刻 状态 量 的 预 估 值 ; #1 为 41 时 
刻 状态 量 的 最 优 估计 值 ; x, 为 点 历 元 时 刻 状态 量 的 一 阶 导数 ; At 为 相 邻 历 元 之 间 的 时 间 间 隔 ; x, 为 
t, 历 元 时 刻 状 态 量 的 二 阶 导 数 ， 
2.2.3 求解 状态 量 的 最 优 估 计 值 x， 
经 上 述 步 又 ， 已 经 求 得 状态 量 值 和 状态 量 的 预 估 值 ,总 历 元 时 刻 直 接 求 解 获 得 的 状态 量 x, FIL, 历 
元 时 刻 获 得 的 状态 量 的 预 估 值 过 ， 现 在 需要 把 它们 组 合 起 来 ， 求 得 二 历 元 时 刻 状态 量 的 最 优 佑 计 值 %， 即 
Z, =a, + a4t, + a4, (6) 
Reb, avs ass a; 可 以 由 广义 融合 最 优化 求解 模型 ( (3) 式 ) 采 用 直接 法 等 优化 算法 求解 得 到 。 将 得 到 
的 w a, a, 代入 (6) 式 可 得 到 4 历 元 时 刻 状态 量 的 最 优 估计 值 & 。 
2.2.4 求解 一 组 优化 的 最 优 佑 计 值 解 
对 上 述 4 个 步骤 不 断 进行 迭代 、 递 推 , 直 到 结束 ， 就 能 够 得 到 一 组 优化 的 最 优 估计 值 解 x,， 
2. 3... kdo 
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3 广义 融合 模型 之 中 权 值 分 析 


在 广义 融合 方法 数学 模型 ((3) 式 ) 中 a, + ast, t au 7 xL RRE BRE GS CRT B 925 ; 
a, + ast, + as) — x, 表示 逼近 函数 与 状态 量 的 预 估 值 之 间 的 残 差 。 它 们 的 权 系数 k, k 的 取 值 大 小 会 
影响 它们 各 上 自 的 残 差 对 目标 函数 的 贡献 ， 所 以 权 值 的 选择 极为 重要 。 故 在 确定 权 系数 之 前 要 分 析 状 态 
量 的 预 估 值 x, 、 状 态 量 的 测量 值 x, 与 逼近 函数 之 间 的 残 差 大 小 比例 关系 ; 在 融合 数据 时 ， 要 考虑 被 
融合 数据 自身 的 精度 及 误差 特性 等 诸多 因素 。 通 过 改变 、h 的 大 小 可 以 调控 广义 融合 最 优 估计 值 
与 状态 量 及 预 佑 值 的 靠近 程度 。 

用 仿真 验证 方法 的 可 行 性 ， 即 生成 两 组 数据 : 一 组 数据 随机 误差 较 大 但 偏 移 较 小 (随机 误差 
-3-3， 偏 差 为 1) ; 为 一 组 随机 误差 较 小 但 有 和 较 大 的 偏 移 ( 随 机 误差 -1-1， 仿 差 为 3)。 数 据 生成 后 代 
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入 广义 融合 求解 模型 ， 选 取 不 同 的 4 对 权 值 进行 组 合 分 析 ， 所 得 结果 如 图 2。 
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图 2 (a) k 21, k, 20 时 最 优 估计 值 偏 移 ; (b) k=0.7、k,=0.3 时 最 优 估 计 值 偏 移 ; (c) k 20.3, 
k,-0.7 时 最 优 估计 值 偏 移 ; (d) k, 20, k, 21 时 最 优 估 计 值 偏 移 
Fig.2 (a) k, 21, k,=0, offset of optimal estimation; (b) k=0.7, k,=0.3, offset of optimal estimation; 


(c) k 20.3, k, 20. 7, offset of optimal estimation; (d) k; 0, k,=1, offset of optimal estimation 


从 图 2(a) (b), (e). (d) 4 幅 图 可 以 看 出 : 当 逐渐 减 小 而 h 逐渐 增 大 时 ， 状 态 量 的 最 优 估 
计 值 开始 远离 状态 测量 值 向 预 估 值 交 近 ; 随 着 越 来 越 小 ，%, 越 来 越 大 ， 这 种 现象 更 加 明显 。 在 数 
据 处 理 过 程 中 ， 如 果 有 一 组 数据 序列 的 某 些 测量 量 发 生 突 跳 导致 数据 不 准确 时 ， 便 可 以 通过 改变 权 值 
组 合 来 降低 这 些 突 跳 数 据 对 获取 最 优 估计 值 的 贡献 ， 使 获得 的 最 优 估 计 值 更 加 平稳 。 


4 广义 融合 方法 在 卫星 导航 中 的 应 用 及 仿真 与 实测 分 析 


4.1 广义 融合 在 卫星 导航 中 的 应 用 

以 卫星 定位 导航 为 例 说 明 采 用 广义 融合 方法 的 实用 性 : 卫星 导航 定位 解 算 出 用 户 的 位 置 坐标 为 
x. y. z 3 个 方向 ， 本 文 只 对 一 个 方向 的 用 户 位 置 作 详 细 分 析 ， 其 余 两 个 方向 使 用 相同 方法 修正 即 可 
获得 用 户 位 置 的 完整 坐标 。 
4.1.1 获得 六 时刻 位 置 最 优 估 计 值 解 x， 和 相应 时 刻 的 速度 

将 从 初始 时 刻 到 到 历 元 时 刻 接 收 机 的 位 置 量 设 为 x,(i=1, 2,…, m)。 通 过 对 历 元 时 刻 以 前 的 
位 置 测量 量 采 用 分 段 最 小 二 乘 拟 合 的 方法 (5 个 历 元 时 刻 点 做 一 次 拟 合 ) 对 位 置 测 量 量 进行 处 理 ， 用 处 
理 后 通过 通 近 多 项 式 (6) 获 得 的 上 时 刻 的 值 ， 作 为 上 时 刻 位 置 的 最 优 估计 值 解 x., 52 。 由 于 在 导航 定 
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位 中 速度 的 测量 值 精 确 度 较 高 ， 初 始 时 刻 的 最 优 速度 值 及 之 后 各 时 刻 的 最 优 速度 值 都 直接 选用 对 应 时 
刻 的 速度 测量 值 即 可 。 
4.1.2 获取 14, 时 刻 位 置 量 的 预 估 值 x,,， 

由 第 1 步 已 经 得 到 了 时刻 位 置 的 最 优 估 计 值 解 x, 和 相应 时 刻 速度 的 最 优 值 v,， 利 用 (5) 式 
解 算出 Li 时刻 位 置 量 的 预 估 值 %,,,。(5) 式 中 %%, 为 上 历 元 时 刻 状态 量 的 一 阶 导数 即 ,时 刻 的 速度 
值 ; x, X t, 历 元 时 刻 状态 量 的 二 阶 导 数 即 加 速度 值 ， 加 速度 本 文 暂 不 考虑 。 

4.1.3 获得 上 ,时刻 最 优 估计 值 解 x ， 

通过 上 述 步 又 ,已 经 得 到 包 , 时 刻 位 置 量 的 预 估 值 ， 相 应 时 刻 位 置 量 的 测量 值 可 以 通过 测量 数据 
直接 获得 。 因 为 实际 的 测量 值 会 因为 误差 等 因素 导致 最 终结 果 不 精 确 ， 所 以 本 文 在 测量 值 与 预 估 值 之 
间 添 加 一 个 浆 值 。 当 两 者 之 间 的 绝对 差 值 超过 这 个 浆 值 时 ， 便 认为 测量 值 已 经 发 生 了 跳 变 ， 此 时 (3 ) 
式 中 的 权 值 组 合 二 0、h, 二 1， 使 最 优 估 计 值 的 解 向 预 估 值 方向 通 近 ;， 当 两 者 之 间 绝 对 差 值 没有 超过 
闵 值 时 ， 认 为 测量 数据 在 误差 范围 之 内 ， 就 要 选择 合适 的 权 值 组 合 ， 再 代入 (3) 式 ,组 合 求解 获得 最 
优 估计 值 xi。 

4.1.4 获得 经 优化 的 一 组 最 优 估 计 值 解 
上 述 第 2、 第 3 步 不 断 闪 代 、 递 推 、 融 合 就 能 获得 一 组 优化 的 最 优 估计 值 解 忆 4 xus. n. Kondo 


42 f 
lis 4.2.1. 广义 融合 对 突 跳 数据 的 修正 = : 
— 卫星 定位 系统 的 精度 主要 取决 于 两 方面 : 一 «l| 一 e 一 广义 融合 最 优 估计 值 | 
e 是 观测 卫星 在 空间 的 分 布 情况 ， 通 常 称 为 卫星 分 
一 布 几 何 衰减 因子 ; 二 是 各 观测 量 的 精度 00 。 大 s 
FK 气 层 干扰 、 对 流 层 延 迟 改正 后 的 残 差 、 星 历 误差 。 。? 
一 及 多 路 径 效应 等 影响 导致 最 后 的 测量 数据 产生 偏 > 0| oe, a pe 
c2 差 或 是 突 跳 ， 而 广义 融合 方法 能 够 很 好 地 修正 这 
些 误差 数据 。 现 生成 一 组 随机 误差 在 (一 0.5- 
0.5) 的 数据 ， 分 别 在 第 10、20、30、40 4 个 点 4 
处 加 上 较 大 的 偏 移 量 构成 突 跳 点 ， 然 后 用 广义 融 à 
合 方法 对 这 组 数据 进行 处 理 ， 仿 真 结果 如 图 3。 " 2 T B 


t- 历 元 /s 


由 图 3 可 以 看 出 ， 广 义 融合 法 对 4 个 突 跳 点 
的 数据 都 做 了 实时 修正 ， 其 结果 不 但 在 没有 发 生 
突 跳 的 部 分 很 好 地 遵循 了 原始 测量 值 的 变化 规律 ， 
而 在 有 突 跳 的 时 段 ， 也 能 有 效 排除 突 跳 的 影响 。 
4.2.2 最 小 二 乘 拟 合法 与 广义 融合 法 对 误差 数据 处 理 对 比 

在 卫星 数目 多 于 4 颗 ， 各 测量 值 之 间 相 互 独立 且 符 合 高 斯 正 态 分 布 的 情况 下 ， 采 用 最 小 二 乘 算法 
解 算 接收 机 的 三 维 坐标 是 最 经 典 的 办 法 。 因 为 最 小 二 乘 算 法 能 使 各 观测 伪 距 之 间 的 残 差 平方 和 最 小 ， 
从 而 使 获得 的 估计 值 解 达 到 最 优 。 但 是 多 路 径 效 应 等 因素 产生 的 误差 却 时 刻 影响 其 定位 精度 ， 这 导致 
各 实际 测量 值 之 间 并 不 是 相互 独立 ， 甚 至 有 些 情 况 下 不 符合 正 态 分 布 ， 这 时 若 仍 采用 最 小 二 乘 算法 就 
很 难保 证 最 后 能 获得 最 优 估计 值 解 ' 站 ， 而 广义 融合 法 在 此 种 情况 下 依然 适用 。 

下 面 通过 仿真 实验 验证 此 方法 的 有 效 性 。 令 测量 值 的 理论 值 为 50， 这 时 在 理论 值 的 基础 上 加 串 
声 生 成 一 组 带 有 随机 误差 的 数据 ， 作 为 测量 仿真 计算 数据 ， 为 了 便于 体现 广义 融合 法 的 优越 性 ， 分 别 
用 广义 融合 法 和 最 小 二 乘 算法 对 这 组 数据 进行 处 理 ， 修 正 后 的 效果 如 图 4。 

经 计算 ,得 到 最 小 二 乘 拟 合 值 的 标准 差 为 oa=3. 626 396 349 003 966， 广 义 融 合 最 优 估 计 值 的 标准 
差 为 ao=0.758 068 193045 756， 广 义 融 合法 稳定 程度 要 比 最 小 二 乘 拟 合法 提高 大 约 79%( 图 4); 而 且 
广义 融合 法 的 偏差 也 较 小 ， 误 差 数 据 的 修正 值 更 加 靠近 理论 值 。 


图 3 广义 融合 对 突 跳 数 据 的 修正 


Fig.3 Correction of the data of sudden jump 


with generalized fusion 
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4.2.3 实测 数据 分 析 


(1)2015 4E 5 H 11 日 在 中 国 科学 院 国家 天 文 台 用 中 科 微 ATGM332D 卫星 定位 芯片 的 接收 机 做 动 


态 测 试 。 从 最 终 获取 的 接收 机 三 维 坐 标 中 截取 了 方向 上 一 
二 乘 拟 合法 对 其 进行 处 理 ， 结 果 如 图 5。 


刻 的 原始 速度 值 ， 分 别 用 广义 融合 法 和 最 小 
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图 4 最 小 二 来 拟 合法 与 广义 融合 法 对 
误差 数据 处 理 的 对 比 
Fig.4 Comparison of error data processed by the 
least squares fitting method and the 


generalized fusion method 


由 图 5 可 以 看 出 最 小 


Y/m 


段 带 有 明显 突 跳 的 原始 定位 数据 和 对 应 时 


测量 值 
—e— 广义 融合 最 优 估计 值 
一 一 最 小 二 乘 拟 合 值 


图 5 广义 融合 法 和 最 小 二 乘 拟 合法 对 
实测 数据 处 理 的 对 比 
Fig. 5 Comparison of measured data processed by 
both the generalized fusion method and the 


least square fitting method 


二 乘 拟 合法 对 某 些 突 跳 点 并 没有 很 好 的 修复 作用 ， 只 是 在 少数 点 处 遵循 了 原 


始 数据 的 变化 规律 ， 但 是 整体 数据 偏离 了 实测 数据 ; 而 广义 融合 法 能 实时 修正 原始 数据 中 的 坏 数据 ， 
并 且 在 没有 发 生 突 跳 的 部 分 很 好 地 跟随 原始 数据 的 运动 轨迹 。 

(2) 当代 社会 对 时 频 基准 的 要 求 越 来 越 高 ， 氢 原子 钟 、 饮 原子 钟 等 多 种 高 精度 时 频 基准 已 被 相继 
研发 出 来 。 然 而 还 可 以 通过 另外 一 种 方法 获得 更 高 精度 和 稳定 度 的 时 频 基 准 ， 那 就 是 组 合 处 理 和 精细 


加 工 原 子 钟 组 的 输出 数据 ， 获 得 性 


能 更 优 恨 的 时 频 输 出 数据 。 本 文 把 氧 原子 钟 与 馅 原子 钟 的 输出 数据 


用 广义 融合 法 融合 起 来 。 图 6(a) 显示 了 饮 原 子 钟 与 氨 原 子 钟 输出 的 频率 数据 变化 图 ， 通 过 对 比 可 以 


发 现 饮 原 子 钟 频率 的 随机 误差 偏 大 ， 但 长 期 稳定 度 很 好 ; 而 氢 原 子 钟 则 不 同 ， 
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它 随 机 误差 的 幅度 很 小 ， 
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图 6 (a) 氢 钟 、 铬 钟 时 频数 据 对 比 ; (b) 广义 融合 前 后 时 频数 据 对 比 
Fig.6 (a) Frequency data of hydrogen and cesium clock; (b) Comparison of time 


frequency before and after the generalized fusion 


c— 


34 天 文 研究 与 技术 13 郑 


但 是 随 着 时 间 的 增长 其 频率 出 现 很 明显 的 漂移 。 广 义 融 合法 能 把 这 两 种 原子 钟 的 频率 特性 的 优势 组 
合 起 来 ， 从 而 获得 随机 误差 小 并 且 长 期 稳定 度 好 的 时 频 基 准 。 用 这 一 方法 对 国家 授时 中 心 的 饱 原子 
钟 和 氧 原子 钟 输出 的 时 频数 据 进 行 融合 处 理 ， 用 氢 钟 时 频 的 变化 率 修 正 饮 钟 的 时 频 输 出 ， 处 理 结果 
见 图 6(b) 。 


5 结 it 

本 文 针对 卫星 导航 定位 数据 中 出 现 的 坏 数据 ， 尤 其 是 突 跳 数 据 的 修正 ， 在 广义 延 拓 插值 外 推 法 的 
基础 上 提出 了 一 种 新 的 数据 融合 方法 一 一 广义 融合 法 。 通 过 仿真 和 实测 数据 分 析 ， 与 最 小 二 乘 算法 进 
行 对 比 后 得 出 : 广义 融合 法 较 最 小 二 乘 拟 合法 更 能 够 修复 定位 数据 中 产生 的 突 跳 数据 ;广义 融合 法 比 
最 小 二 乘 拟 合法 对 数据 的 处 理 可 靠 性 更 高 ;广义 融合 法 实用 性 好 ， 数 据 平稳 性 好 ， 外 推 稳定 度 高 ， 有 
望 成 为 一 种 新 的 解决 数据 融合 问题 的 普 适 求解 方法 。 
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A New Method of Data Fusion—— the Generalized Fusion Method 
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Abstract: The purpose of data fusion is to merge or synthesize two data sets so as to make them become 
a better and complete data set. Based on extended inter-polation method, generalized fusion method is a 
method which strives to create a mathematical model that consists of two or more sets of equation and data sets. 
This method is not only able to bring out its full potential of data and measurement resources, but also capable 
of using innovative integration and mutual algorithm to solve the data fasion problem. After simulation and real 
measurement data analysis, conclusion can be drawn that compared with the least square algorithm, the 
generalized fusion method is more suitable for the repair of the data generated in the process of data fusion. 
This paper also introduces the main steps of the establishment and the solution of the generalized mathematical 
model, analyzes the distribution principle of weight in the model and its application in the field of satellite 
navigation data processing. Experimental results show that the method is practical and stable, and it is 
expected to be a new method for solving the problem of data fusion. 
Key words: Data fusion; Generalized fusion; Generalized extrapolation; Model solving; Weight distribution; 


Satellite navigation data processing 


